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For a mixture of alkali-earth atomic gas in the long-lived excited state 3 Po and the ground state 
1 So, in addition to nuclear spin, another “orbital” index is introduced to distinguish these two inter¬ 
nal states. In this letter we propose a mechanism to induce Feshbach resonance between two atoms 
with different orbital and nuclear spin quantum numbers. Two essential ingredients are inter-orbital 
spin-exchanging process and orbital dependence of the Lande g-factors. Here the orbital degrees 
of freedom plays similar role as electron spin degree of freedom in magnetic Feshbach resonance 
in alkali-metal atoms. This resonance is particularly accessible for 173 Yb system. The BCS-BEC 
crossover in this system requires two fermion pairing order parameters, and displays significant 
difference comparing to that in alkali-metal system. 


Magnetic Feshbach resonance (MFR) is a powerful tool 
to tune interaction between atoms to strongly interact¬ 
ing regime, and plays a crucial role in cold atom physics 
Q.J. For instance, for alkali atoms, each atom has an elec¬ 
tronic spin S = 1/2. The interaction between two atoms 
have different potentials depending on whether the total 
electronic spin is singlet or triplet. Thus, one can utilize 
the Zeeman energy to control their relative energy and to 
reach a scattering resonance. However, alkali-earth atom 
(like Sr) or alkali-earth-like atom (like Yb) has fully oc¬ 
cupied outer shell, and their total electron spin is zero. 
Thus, it is conventional wisdom that there is no MFR in 
alkali-earth atoms at ground state. Instead, one can tune 
interaction by optical Feshbach resonance j2], but such a 
scheme suffers from strong atomic loss and heating. 

Another significant feature of alkali-earth atom is the 
existence of a long-lived excited state 3 Pen in which one 
electron is excited to p-orbital and the total electronic 
spin S = 1. The dipole transition to ground state 1 Sq 
is “spin-forbidden” and therefore the lifetime of excited 
state can be as long as a few seconds. This is used for 
atomic clock transition. Considering an atomic gas mix¬ 
ture of 3 P 0 (denoted by |e)) and 1 S 0 (denoted by |g)) 
states, in addition to nuclear spin degree of freedom, one 
introduces another so-called “orbital” degree of freedom 
to label the internal state of atoms GHH]. Previously at¬ 
tentions have been paid to MFR between 3 P 2 state and 
1 So state due to anisotropic interactions, but these res¬ 
onances are generally quite narrow [6]. Moreover, since 
J = 0 for 3 P 0 , even such a MFR does not exist between 
3 P 0 and 1 Sq. 

In this letter we propose an alternative mechanism that 
can lead to a Feshbach resonance (FR) between atoms in 
|e) and \g). Though this FR is also controlled by mag¬ 
netic field, the mechanism of how it works is quite differ¬ 
ent from MFR in alkali atoms. Two essential ingredients 
are inter-orbital (nuclear-)spin-exchanging interactions, 
which has been observed in recent experiments EHS], and 
the small difference in the nuclear Lande g-factor Sg be¬ 



FIG. 1: (a) Original energy level diagram. A E denotes the 
excitation energy between |e) and | g). 8 g = Bg g p,B and 
8 e = Bg e /iB are Zeeman energies of | g) and |e) states, re¬ 
spectively. Two states in open channel are occupied. Arrows 
indicate an inter-orbital (nuclear-) spin exchanging scatter¬ 
ing process, which couples open channel |g f; e i") and closed 
channel \g ■Re f). S = 5 e — 5 S = B(5g)g.B is the Zeeman 
energy difference between two channels and 8g = g e — g g is 
the difference in Lande g-factor. (b) Reorganized energy level 
diagram for many-body Hamiltonian, in which open channel 
states appear in lower energy. 


tween different orbital states (|e) and |.g)) [7]. The orbital 
degree of freedom plays the role as electronic spin in FR 
of alkali atoms, and we therefore name it as “Orbital 
Feshbach resonance” (OFR). 

Two-body Problem with Pseudo-Potential Approach. 
For simplicity, we first illustrate the basic idea with 
pesudo-potential. And without loss of generality, we take 
two nuclear spin states (mj and mj + 1) denoted by | t) 
and | |). Their energy level diagram is shown in Fig. 
[lja) and explained in the figure caption. We consider 
two atoms, either one in \g j.) and the other in |e f) (de¬ 
noted by open channel |o) = (| g e f) — |e t; 5 4-))/v / 2), 
or either one in |g,f) and the other in |e |) (denoted by 
closed channel |c) = (|g f> e I) — |e i;g t))/ a/ 2)- The 
threshold energy of these two channels differs by S due 
to the difference in the Lande g-factor. For the rela- 
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tive motion between these two-atoms (with mass m), the 
non-interacting Hamiltonian is written as 

/ fi 2 V 2 \ h 2 V 2 

H 0 =( -— + 6) |c)<c|-—|o)(o|. (1) 

\ m J m 

The interaction part depends on whether orbital degree 
of freedom forms a singlet or triplet. For s-wave scatter¬ 
ing we introduce two anti-symmetrized bases 

l±> = \{\9e) ± Ie#))(| n) T I It)) = ^=(\c) T |o», (2) 

the Huang-Yang pseudo-potential is diagonal in this 
bases with two different scattering lengths and a ~, 



When rotated into the |o) and |c) bases, the interaction 
potential V becomes 


V = Vo (|o)(o| + |c)(c|) + Vi (|c)(o| + |o)(c|), (4) 

where Vj = and a s o denotes (a+ + 

a“)/2, and a s i denotes (a“ — a+)/2. The V) term de¬ 
scribes an inter-orbital spin exchanging process, as illus¬ 
trated in Fig. [I] which couples the open and closed chan¬ 
nels. A positive a s o is always associated with a bound 
state with binding energy £& = —h 2 /(ma 2 0 ). Therefore, 
when S ~ £&, one will expect a scattering resonance in 
the open channel. 

The two-body wave function can be written as 


mS / h 2 —k 2 r 

\o)+m- —-—ic). 

A A (5) 

Solving the Schrodinger equation (H 0 + V)ip = Eip with 
E = h 2 k 2 /m , one can find 

I mS 

(1 + ika s0 )f o (k) - a s i y -p- - k 2 f c {k) + a s0 = 0; (6) 

ika sl f 0 (k) + ^1 - a s0 \j~^ - k^j fc(k) + a s i = 0. (7) 


i> = 


„ikr 


+ fo{k) - 

r 



FIG. 2: Scattering length a s between |e f) and |g 4-) as a 
function of magnetic field for 1,3 Yb atom. Blue dashed line 
and red dash-dotted line is the real and imaginary part of a s 
obtained from the zero-range pseudo-potential. Solid line and 
red dotted line are a s obtained from the finite range potential. 
Here we take ( Sg)fj,B = 2irhx 112Hz/G, a| = 3300ao— *0.78ao, 
&S = 219.5ao [U®[8; and ro is taken as van der Waal length, 
which equals to 84.8ao 0, with a o the Bohr’s radius. 


For r < ro, the Hamiltonian is written as 
/ h 2 \7 2 \ 

H= E + W(r) |t)<i|, (10) 

«=+,- » ' 


where we have assumed that 5 is much smaller than en¬ 
ergy scale of short-range potential V l (r) such that it can 
be safely ignored in this regime. Each V l (r) (i = +, —) 
corresponds to an s-wave scattering length a*, that is to 
say, the wave function '0* = u l [r)/r in r < ro regime sat¬ 
isfies the boundary condition u l, (r)/u l (r) | r=T . 0 = l/(ro — 
a l s ). In r < ro regime the wave function can be written 
in a general form 


u(r) = M + (r)|+) + Au (r)|—) 

i+. 

Ic) + 


u 1 (r) + Au (r). . — u + (r) + Au (r) 


V2 


V2 


( 11 ) 


Straightforward calculation yields the scattering length 
a s in the open channel as 


a s = ~fo{k = 0) 


-q s0 + \/mS/H 2 (a 2 0 - a 2 ^ 

a s o\J mS / h 2 — 1 


Two-body Problem with Finite Range Potential. We 
can further more rigorously demonstrate the OFR with 
a coupled two-channel model with finite range ro- When 
r > ro, two atoms are non-interacting, and the zero- 
energy s-wave wave function ip = u{r)/r with 

j |c) +/3(r - a s )|o). (9) 



By matching boundary conditions between wave func¬ 
tions Eq. [9] and Eq. 11 at r = ro for |o) and |c) channels 
independently, and utilizing the boundary condition for 
each u l , one can obtain A and a s , where 


-a s0 + yJmfjW [(a 2 0 - o^) - r 0 a s0 ] 


y/mS/h 2 (a s0 - r 0 ) - 1 


( 12 ) 


In the limit ro —> 0, Eq. 12 recovers the result of Eq. [8] 
Thus we have demonstrated the OFR phenomena. Eq. 
and Eq. 


[8] and Eq. 12 also show that the difference between a” 
and a+, i.e. a s i 0, is crucial, as a s becomes a constant 
as a s o if one sets a s i = 0 in Eq. [8] and Eq. 12 


u(r ) = a exp 
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FIG. 3: Pairing order parameters (a) and chemical potential 
(b) as a function of 5/5 Tea across an OFR, with 5 lea given 
by Eq. [13] (ro = 0). In (a), red solid and blue dashed lines 
are pairing order parameters in the open and closed chan¬ 
nel, respectively. For different 5/5 rea , corresponding values of 
1/(A:fo s ) are marked by arrows, where a a is the open chan¬ 
nel scattering length. Here at and a 7 are taken from 1,3 Yb 
atom, as described in caption of Fig. [5J and n = fcp/(37r 2 ) 
is taken as 5 x 10 13 cm~ 3 in our calculation. For comparison, 
black dotted lines in (a) and (b) show pairing order parameter 
and chemical potential for single-channel BCS-BEC crossover 
with same l/(IcFOa). 


Orbital Feshbach Resonance. Eq. [12] shows that a s 
diverges when 


S = S les 


ti 2 

m(a s0 - r 0 ) 2 ' 


(13) 


This determines the position of OFR. It also indicates 
that the precise location of OFR may be altered by short- 
range details. In addition, we also obtain that a s has a 
zero-crossing when 


<5 = 5 0 = 


h 2 

m(a s0 - r 0 - a^/aso) 2 ' 


(14) 


Since S = B(5g)g.B, we take at, a~, ro and (Sg)g .b from 
173 Yb atom measurement reported so far, and plot a s as 
a function of magnetic field B , as shown in Fig. i 173 Yb 
is a unique system to observe an OFR, since a 3- is as large 
as ~ 10 3 ao- Therefore, it does not require large magnetic 
field to reach the OFR. Otherwise, if the scattering length 
is small, it requires a much high magnetic field (10 3 — 


10 4 G for s7 Sr) to reach this resonance. On the other 
hand, since the essential ingredient of inter-orbit spin¬ 
exchanging process also exists between 3 Pi (or 3 P 2 ) state 
and 1 S 0 , OFR can also exist in these mixtures, where 
the g-factor difference is quite large, and it does not rely 
on a large scattering length at zero-field. Nevertheless, 
the scattering behavior between these states is also more 
complicated, and the details are left for future studies. 

Regarding loss nearby OFR, there are two mechanisms. 
First, the inelastic two-body loss, described by the imag¬ 
inary part of scattering length, is also shown in Fig. [2] 
Secondly, the shallow bound state is also subjected to 
three-body loss, which cannot be rigorously captured by 
our two-body calculation. Nevertheless, we phenomeno¬ 
logically include this effect by introducing a decay term 
on the shallow bound state. We find that it will not 
affect the real part of scattering length, as long as the 
bound state has a reasonable lifetime, but only enhance 
its imaginary part (10l . 

Here we would like to contrast the OFR in alkali-earth 
atom with MFR in alkali-metal atom. In MFR, when 
two atoms interact at short distance, the interaction po¬ 
tential are different for total electronic spin singlet and 
triplet. While for OFR, the interaction potentials are 
distinguished by orbital singlet or triplet. In MFR, the 
coupling between two channels is due to hyperfine inter¬ 
action. While in OFR, orbital dependent Lande g-factor 
can be viewed as coupling between orbital and nuclear 
spin, because of which two channels are coupled. In this 
analogy, orbital degree of freedom in OFR plays the same 
role as electronic spin degree of freedom in MFR. 

Two Order-Parameters BCS-BEC Crossover. Since 
(Sg)fi b is five orders of magnitude smaller than the g e p b, 
comparing OFR with a MFR in the same magnetic field 
regime, the energy separation between open and closed 
channel is much larger in the MFR case than that in the 
OFR case. In a MFR case, this energy separation is a 
few orders of magnitude larger than the Fermi energy. 
Therefore, in a BCS-BEC crossover theory studied be¬ 
fore HU, one can either start with a single-channel model 
only, or with a two-channel model but only including the 
bound state of the closed channel. The scattering states 
in the closed channel are never important. However, the 
situation in the OFR case is considerably different. Con¬ 
sidering a typical density of Fermi gas, S is comparable 
or can be even smaller than the Fermi energy. Thus, we 
have to take into account scattering states in both open 
and closed channels. This requires introducing two self- 
consistent paring order parameters for open and closed 
channels, respectively. Below we shall present such a for¬ 
malism for crossover across an OFR. 

We consider the situation that in the non-interacting 
limit, total N fermions are equally populated in the two 
states in the open channel (|e f) and |g 4-))- We note 
that both N e = N e - j- + N : e j. and Nf = N g -j- + IV e f are good 
quantum numbers. Subtracting the Hamiltonian with a 
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constant term (A E + S/2)N e — (S e + S g )Nf/2, one can 
show the level diagram can be reorganized as shown in 
Fig- 0b), such that the two states in the open channel 
appear in the lower energy, with an energy separation of 
<5/2 below the two states in the closed channel. Then in 
the non-interacting limit (say, S = So), S/2 is larger than 
the Fermi energy so that only open channel is equally 
populated. The many-body Hamiltonian can be written 
as 


H — H 0o + Hoc + -f —— A}_A_, (15) 

Hqo = 5Z £k ( C ff4-k c stk + C e|k c eTk) (16) 

k 

Hoc = ^ £ k + ^ J ( C gtk C flTk + C e4.k C etk), (1") 

where £k = fi 2 k 2 /(2m) — g, and 

A+ = ^ —kCf;^l.k Cg^—kCe^k) (1$) 

k 

A- = y> gt —k^e^-k A ('-(jl kX'f’tk)- (19) 

k 

Now we defined two order parameters as A + = 
g + (A)/2 and A_ = g_(A)/2, and we can perform mean- 
field decoupling of the interaction term which leads to 


Hmf —Hqo + Hoc + (A + A + + h.c.) 


+ (A_i_ + h.c.) 


2|A + | 2 

9+ 


2|A| 2 

9- 


( 20 ) 


Following the standard BCS theory to diagnoalized Hmf 
with the Bogoliubov transformation and minimizing the 
ground state energy with respect to both A + and A_ 
[12] . we reach two coupled gap equations 


i-H 

o| u 

<l< 

<1< 

+ 
i—i 

= E 

i 


2m 

47rft 2 aj" 

4nh 2 a s 


A c 2 

h 2 k 2 

m 

m J 

k 

V ( £k + 2 ) + 


Ac _ i 

A 0 1 

1 _L Ac ' 
Ao 

ii 

1 

2 TO 

(21) 

47 Th 2 a+ 
m 

4nh 2 aJ 

m 

V^ + IAol 2 

h 2 k 2 


where A 0 = A_ — A + and A c = A_ + A + are pairing 
order parameters in the open and closed channels, re¬ 
spectively. Here, in contrast to usual BCS-BEC crossover 
where scattering length is the tunable control parameter, 
here both a/ and a~ are fixed. Instead, S is the tunable 
parameter to control the crossover. Moreover, the num¬ 
ber equation also includes contribution from both chan¬ 
nels 


^ = E 

k 


(0 

£ k + | j 

V V £ t + A 0 2 

\J ( £ k + f) + A c 2 / 


( 22 ) 


Solving gap equation Eq. [21] with Eq. [22] we find A 0 
has a 7r phase difference from A c , and we determine both 
|A 0 |, | A c | and g as a function of <5, as shown in Fig. 0] 
We find (i) when S <5 res , the system is away from OFR. 
In this case we find small |A 0 |/Ef, |A C |/Ep and g —>■ Ep, 
which is the typical behavior in the BCS regime, (ii) As 
<5 —> S res , both | A 0 1 and |A C | increase rapidly toward 
the same order as Bp and meanwhile, g decreases. This 
feature is qualitatively consistent with a crossover from 
BCS to the unitary regime, (iii) While when S < S res and 
<5 —> 0, both pairing gaps saturate instead of continuously 
increasing toward deep BEC limit. This is consistent 
with that a s finally saturates to a s o- 

We also plot the ratio |A 0 /A C | in the inset of Fig. 0a), 
we find when S <5 res , this ratio decreases toward zero, 
and thus the L.H.S. of Eq. [0] diverges as it depends on 
A c /A 0 . Effectively, if one compares Eq. 21 with a single¬ 
channel BCS-BEC gap equation |T2], this is equivalent to 
that an open channel scattering length decreases toward 
zero. When <5 ->■ S res , |A 0 /A C | -5> (a+-a/)/(a+ +a~) (~ 
0.875 for the scattering lengths we use). The L. H. S. of 
Eq. [0] approaches zero, indicating an divergent effective 
scattering length. Finally when <5 —>• 0, two channels 
become degenerate and thus this ratio approaches unity. 
We remark that this mean-field calculation does not use 
the results from two-body calculation above. It is an 
independent many-body calculation, while the results are 
qualitatively consistent with two-body results. 

On the other hand, for the typical density we con¬ 
sider, the quantitative behavior is quite different from 
the single channel BCS-BEC crossover. We perform a 
single-channel mean-field calculation, in which we only 
keep the two states in the open channel and use open 
channel scattering length a s (S) given by two-body result 
of Eq. 0 The result is shown by the dotted line in Fig. 0 
and compared with the two-gap theory presented above. 
Remarkably, we find that in the BCS regime, the pairing 
in the two-gap theory is stronger than that in the single 
channel model. While if we lower the density so that the 
Fermi energy becomes much lower than S/2, the results 
gradually converges to the single channel result [~L8l . 

Outlook. Our predication of OFR opens an avenue 
for studying strongly interacting physics in alkali-earth 
atomic gases. The two-gap Fermi superfluid is reminis¬ 
cent of two-gap superconductor. Further studies includ¬ 
ing Gaussian fluctuations can reveal effect of finite life 
time of Cooper pairs, superfluid transition temperature 
and the Leggett mode (relative phase mode between two 
gaps) in strongly interacting regime, which is left for fu¬ 
ture investigations. Moreover, by coupling |e f) and |g J,) 
states with a laser, one can create spin-orbit coupling 
between them, which avoids heating from spontaneous 
emission as in Raman scheme. Our OFR increases at¬ 
tractive interaction between them and can help to reach 
a topological superfluid phase in this system. 
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